Quantum-orbit analysis of high-order harmonic generation by resonant plasmon field 

enhancement 



T. Shaaran 1 , M. F. Ciappina 1 and M. Lewenstein 1 ' 2 
1 ICFO-Institut de Ciences Fotdniques, 08860 Castelldefels, Barcelona, Spain 
2 ICREA-Institucio Catalana de Recerca i Estudis Avangats, Lluis Companys, 08010 Barcelona, Spain 

(Dated: March 8, 2013) 

We perform a detailed analysis of high-order harmonic generation (HHG) in atoms within the 
strong field approximation (SFA) by considering spatially inhomogeneous monochromatic laser 
fields. We investigate how the individual pairs of quantum orbits contribute to the harmonic spectra. 
We show that in the case of inhomogeneous fields, the electron tunnels with two different canonical 
momenta. One of these momenta leads to a higher cutoff and the other one develops a lower cutoff. 
Furthermore, we demonstrate that the quantum orbits have a very different behavior in comparison 
to the homogeneous field. We also conclude that in the case of the inhomogeneous fields, both odd 
and even harmonics are present in the HHG spectra. Within our model, we show that the HHG 
cutoff extends far beyond the standard semiclassical cutoff in spatially homogeneous fields. Our 
findings are in good agreement both with quantum mechanical and classical models. 



I. INTRODUCTION 

In the context of the interaction of matter with strong 
laser fields, high-order harmonics generation (HHG) pro- 
cess [3, Q has attracted considerable interest, since it 
represents a viable route to the generation of coherent 
radiation in the ultraviolet (UV) to extreme ultraviolet 
(XUV) spectral range with high repetition rate. Due to 
this high repetition rate, as well as high coherence degree 
and wavelength tunability, HHG has found numerous ap- 
plications in various areas of science, such as material 
sciences, life sciences and lithography [3[. 

In addition, HHG has been employed for generating 
ultrashort pulses, and even single attosecond pulses Q. 
This allows even more control over the atomic and molec- 
ular processes. For instant, by superposing the XUV 
pulses on the laser field, one can resolve dynamic pro- 
cesses in atoms or molecules with precision of a few at- 
toseconds [5|. The physics behind the HHG can be un- 
derstood by a simple semi-classical three-step model [6|. 
According to this picture, an electron leaves the atom 
or molecule by tunneling through the potential barrier, 
formed by the atomic potential and the laser electric field, 
to reach the continuum. It subsequently propagates in 
the continuum and driven back by the laser filed toward 
its parent ion or molecule. Finally, upon its return, it 
recombines with the core and leads to the emission of 
energetic photons. 

The threshold intensity for generating high-order har- 
monic in noble gases is above 10 13 W/cm 2 , which is far 
beyond the output power of the current femtosecond os- 
cillators. Nowadays, chirped-pulse amplification (CPA) 
is used to exceed the threshold intensity. CPA is a com- 
plex process and it requires multi-pass amplifier cavities 
in tandem. On the other hand, as far as the applications 
of the HHG are concerned, the optimization of HHG effi- 
ciency and the extension of the cutoff to short- wavelength 
region are important aspects to consider [7|, [8[ . The cut- 
off frequency of the generated harmonics can be extended 
either by reducing (increasing) the laser frequency (wave- 



length) or increasing the peak field amplitude. There is, 
however, limitations in tuning these parameters. In fact, 
reducing the laser frequency leads to a significant drop 
in the harmonic yield due to the longer electron excur- 
sion times 0, [To(| , while increasing the laser field intensity 
produces depletion of the ground state. 

A recent demonstration has shown that the surface 
plasmon resonance could provide a possible solution to 
overcome the problems cited above The local field 
enhancement induced by a resonant plasmon within a 
metallic nanostructure requires no extra cavities or laser 
pumping for amplification. In this scheme, the local elec- 
tric fields can be intensified by more than 20 dB HI, (3, 
an amplification that exceed the threshold laser inten- 
sity for HHG generation in noble gases. In addition, the 
pulse repetition rate remains unchanged without adding 
any additional cavities or extra pumping. Moreover, each 
nanostructure acts as a point-like source to generate har- 
monics radiation, which through constructive interfer- 
ence can provide even more focused coherent radiation. 
This gives a wide range of possibilities to spatially rear- 
range nanostructures to shape or enhance HHG spectral 
and even obtain a new physics. The locally enhanced 
field, has a distinct spatial dependency, which gives an 
enormous extension to the HHG cutoff [14{ . 

HHG based on plasmonics can be understood as fol- 
lows [11]: the external femtosecond low intensity pulse 
couples to the plasmon mode and induces a collective os- 
cillation of free charges within the localized regions of the 
nanostructure. The free charges redistribute the electric 
field around the nanostructure vicinity, in such, to form 
a spot of highly enhanced electric field. The enhanced 
field, which largely depends on the geometrical shape of 
the metallic nanostructure, exceeds the threshold inten- 
sity required for HHG. As a result, by injection of noble 
gases into the spot of the enhanced field, one can gen- 
erate high order harmonics. In Ref. [11[, the output of 
the femtosecond oscillator, which was a pulse with 10 
fs pulse duration, 800 nm of wavelength and intensity of 
10 11 W/cm 2 , was directly focused onto a bow-tie nanoan- 
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tenna. As a result of the laser pulse interaction with the 
nanostructure, the field intensity is enhancement by 2-4 
order of magnitude, which is sufficient enough to produce 
XUV wavelengths from the 7 th (114 nm) to the 21 st (38 
nm) harmonics in Argon. 

Hitherto, the theoretical approaches for studying 
strong field phenomena are largely based on the as- 
sumption that the laser electric field is spatially homo- 
geneous in the region where the electron motions take 
place [HI, [l6j]. This assumption, however, does not hold 
for the field enhanced by resonant plasmons. Indeed, the 
strong confinement of the electrons in the plasmonic hot 
spots generates a spatially inhomogeneous electric field, 
which strongly influences the subsequent motion of the 
electrons in the strong field phenomena. As a result, new 
physics will emerge in the interaction between matter and 
strong laser fields. 

Since the first observation of HHG, different theo- 
retical models, including solving the Time Dependent 
Schrodinger Equation (TDSE) and the Strong Field 
Approximation (SFA), have been applied to describe 
this phenomenon (for details see the review articles 
in [13, [HI). On the other hand, the HHG by resonant 
plasmon field enhancement is a new topic and it has been 
considerably less well studied since it is far more difficult 
to measure and model. Up till now, two experiments have 
been performed to measure HHG based on this kind of 
field 11 [l9| , while only few theoretical approaches have 
been developed [20l-[22| . One should note, however, that 
the interpretation of Ref [ll[ remains controversial [23l — 

HI. 

In our previous paper Q , in which we employed SFA 
to investigate the HHG by resonant plasmon field en- 
hancement, the vector potential field was defined in an 
approximated way starting from the inhomogeneous elec- 
tric field. In the present paper, we improve our model 
even more by applying the same approximation directly 
to the inhomogeneous electric field rather than the po- 
tential. Subsequently, the vector potential is now derived 
from the electric field. In addition, we use a SFA based 
on saddle point methods rather than a full numerical SFA 
approach to obtain the HHG spectra. The saddle point 
methods lead to equations that can be directly related to 
the classical equations of motion of an electron in a laser 
field. As a result, they provide a space-time picture which 
gives us additional physical insight. In this work, we 
scrutinize the individual electron trajectories in compar- 
ison to their classical counterparts and demonstrate their 
contributions to the HHG spectra. In addition, since the 
imaginary part of the saddle point equations can be re- 
lated to the width of the potential barrier through which 
the electron tunnels, we examine the ionization probabil- 
ity of the electron for each trajectory. 

This article is organized as follows. In Sec. II, we 
present the SFA transition amplitude for high-order har- 
monics, starting from common expressions based on ho- 
mogeneous electric fields (Sec. II. A) and, subsequently, 
by showing how we modify it for the case with non- 



homogeneous fields (Sec. II. B). In Sec. Ill, we discuss 
the saddle-point equations and analyze them in terms of 
quantum orbits in parallel to their classical counterpart 
electron trajectories. In the next section, IV, we present 
the HHG spectra based on the analysis given in Sec. III. 
Finally, in Sec. V, we summarize the paper and state our 
main conclusions. 



II. THEORY 

A. Transition Amplitude 

Generally, there are two main assumptions behind the 
strong field approximation (SFA), namely (i) the influ- 
ence of the laser field is neglected when the electrons are 
bound to their target atoms and (ii) the binding ionic po- 
tential is neglected when the electrons are in the contin- 
uum. As a result, the free electrons in the continuum are 
described by field-dressed plane waves, which are known 
as Volkov states [H, [27j. 

1. Homogeneous fields 

In the Lewenstein model [28j, a well established SFA- 
based method to model HHG, it is assumed that the laser 
electric field does not change with respect to the position 
in the region where the electron motion takes place. In 
this spatially homogeneous field, the SFA transition am- 
plitude for HHG reads (in atomic units) 

/OO f't r 

dt dt 1 d 3 kd* rec (k(t))d ion (k(t')) 
-OO J —OG J 

e -iS(n,k,t,t') +c . c . (1) 
with the action 

So(fi,M,0 = / t k + A ( r ^ dr + I p (t -t)-nt (2) 
Jt' ^ 

and the prefactors 

^ on (k(t / )) = (k(t , )l^(t / )l^o) (3) 

d rec (k(t)) = (Ht)\Odi P -eM- (4) 

Thereby, /c, I pi Q , Hi nt (t'), Odi P and e x denote the drift 
momentum of the electron in the continuum, the ioniza- 
tion potential of the of the field- free bound state |0o), 
the harmonic frequency, the interaction of the system 
with the laser field, the dipole operator and the laser po- 
larization vector, respectively. The vector potential A(t) 
of the laser electric field E(t) is defined by 

A(t) = - f E(t')dt'. (5) 
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Physically, Eq. (pQ) describes a process in which an elec- 
tron, initially in a bound state |0o) with energy J p , in- 
teracts with the laser field by H int (t') at the time t' 
and tunnels into a Volkov state |k(t)). Subsequently, 
from time t' to t, it propagates in the continuum and 
is driven back by laser field to its parent ion. At the 
time t, upon its return, this electron recombines with 
the core and emits high-harmonic radiation of frequency 
Q. The second term in Eq. (pQ), which corresponds to 
the continuum-continuum transitions, can be ignored, 
since it contributes very insignificant to the transition 
amplitude of HHG (for details see [29| ). The main 
drawback of SFA is that it is not gauge invariant. As 
a result, the matrix elements describing the ionization 
and recombination, i.e. Eqs (j3j) and (j4|, have different 
form in the length and velocity gauges. This comes from 
the fact, that both the interaction Hamiltonian Hi nt (t') 
and the Volkov wave function |k(t)) are not translation- 
ally invariant. The interaction Hamiltonian is given by 
H\ nt {t) = r • E(f ) and H? nt (t) = [k + A(f )]/2 in the 
length and velocity gauges, respectively. For the Volkov 
wave function, k(t) = k + A(t) in the length gauge and 
k(t) = k in the velocity gauge. In this paper, we work in 
the length gauge and we assume that the electric field is 
linearly polarized along x-axis. Furthermore, we consider 
an hydrogenic Is state for the field- free bound state \4>o)- 
As a result, the Eqs. (|3)) and (g]) yield 



dion(k(£')) °c 



(k(^) 2 + a 2 ) 3 



E{t f ) 



(k(i)) oc 



k(t)x 



(k(t) 2 + a 2 ) 3 



(6) 



(7) 



2. Nonhomogeneous fields 

We will now consider a case in which that assumption 
made in Sec. Ill A II is not any more valid and the laser 
field has a spatially inhomogeneous character, when the 
HHG process takes place. Before discussing the nonho- 
mogeneous case, we examine how the action of the SFA 
is connected to classical electron trajectories for the ho- 
mogeneous field. The laser potential Vl due to the laser 
field E(t) is defined as 



V L =xE(t), 



(8) 



and the Newton equation of motion for an electron in 
this field is given by 



x(t) 



-V X V L 



(9) 



In here, the force is equal to (minus) the laser electric 
field, i.e. x(t) — —E(t). In the SFA, the action is defined 
in terms of the vector potential field A(t) given by (j5]), 
which is the counter part of the velocity x(t). 



For the inhomogeneous case the electric field now has 
the form E(t, x) and the laser potential is Vl = xE(t, x). 
Thus, the Newton equation of motions become 



x(t) = -xV x E{t, x) - E(t, x). 



(10) 



From Eq. ([TP]) , it is clear that x(i) ^ -E(t,x). There- 
fore, the potential field A(t) of the SFA action should 
correspond to the integration of x(t) with respect to t. 

If the spatial dependence of the laser electric field is 
perturbative and linear with respect to position, then 
the field can be approximated as 



E(t,x) ~ E(i)(l + ex), 



(11) 



where e <C 1 is a parameter that characterize the strength 
of the inhomogeneity. 

Indeed, the above approximation corresponds to the 
first term of the actual field of a plasmonic nanostructure 
with spherical shape [30| . By substituting (JTTJ into ([TO]) , 
we have 



x(t) = -E(i)(l + 2ex(i)). 



(12) 



This is the effective laser electric field that the electron 
feels along the trajectory x(t), which describes its motion 
in the continuum. We will call it the electron trajectory 
effective electric field. 

Classically, the electron trajectory can be found by 
solving Eq ([T^]) . In here, we solve it by applying the 
Picard iteration [31( method and restrict ourselves to the 
first order (for more details see [14[ ). Based on the con- 
dition that the electron starts its movement at the origin 
with zero velocity, i.e. x(0) = and x(0) — v(0) — , we 
obtain 

x(t) = 0(t) - p(t ) - A(t )(t - t ). (13) 

with P(t) = f*dt'A(t'). In addition we assume that at 
time to the potential field is zero; thus 



s(t) = J* dt'A(t'). 



(14) 



By using Eqs. flSJ, (fT3)) and ([T2]). the effective vector 
potential along the electron trajectory A tr (t) reads 



A tr (t) = A(t) + 2eA c (t), 



(15) 



where 



A c (t) = J' dt"A(t") - J' dt"A 2 (t"). (16) 

The next step is to modify the general expression of the 
transition amplitude (pQ) for the HHG, in order to include 
the above defined inhomogeneous field. Consequently, we 
have to replace the electric field and vector potential by 
Eqs. ([TT]) and (fT5]h respectively. As a result, the modified 
action yields 

5(0, k,M') = 5'o(^,k,t,t / ) + 2e^A c (r)[k + A(r)]^r 
+2e 2 \\{r)dr (17) 
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where 5q(0, k, t, t') is defined in Eq. (j2j). 



B. Saddle-point equations 

The transition amplitude (pQ) can be computed either 
numerically or using the saddle-point method [32|, [33| . In 
here, we employ the latter procedure since the solutions 
of the saddle point equations are directly related to the 
classical trajectories. Thus, it allows us to investigate 
the quantum orbits in comparison to the classical tra- 
jectories as well as demonstrating their contributions to 
the cutoff and yield of the HHG. This method requires 
obtaining the saddle points where the action (fT7|) is sta- 
tionary, i.e. for which d t S(Q, k, t, t') = dfS(Q, k, t, t r ) = 
dkSiyt^t^t') = 0. In this paper, we use a specified 
steepest descent method called uniform approximation 
to take care of those saddles points which are not well 
separated (for a detailed discussion see Ref. [34J ). The 
stationary conditions upon t, t' and k lead to the saddle- 
point equations 



limit, in which the solutions of (fT8|) are real. Neverthe- 
less, in here, e is a very small parameter and the electron 
will most likely reach the continuum with tunnel ioniza- 
tion. Thus, the solutions of these saddle point equations 
are still expected to be complex. In addition, the maxi- 
mum kinetic energy that the electron gains in the contin- 
uum is not any more 3.17/7 P , where U p = Eq/ (4cj 2 ) is the 
ponderomotive energy. In fact, it depends on the nonho- 
mogeneous character of the field, i.e. of e and A c (t). For 
positive A C1 the electron gains energy, depending on the 
value of the e, larger than 3.17/7 P , while for the negative 
A C: it would be below the conventional value. 

We now examine the drift momentum k of the electron 
at the time of the tunneling Eq. ([T8]) . For that, we 
consider the limit I p 0, where the electron reaches 
the continuum with zero kinematical momentum. As a 
result, Eq. ([T8]) yields to 



k = -i(0 + eA c (0(lTV3) 



(23) 



[k + A(t')] -2eA(t') = -2I P 



J dr[k + A(r)] + 2er](r) = 



(18) 



(19) 



Unlike the homogeneous case, where k = —A(t'), in here, 
k has two different solutions with one exceeding and the 
other lowering the homogeneous drift momenta. The 
strength of the inhomogeneity e and the shape of the 
A c (t') are the responsible of this deviation. 



and 



with 



and 



n=Pi±^C + J p + 2eA(t) 



A(t) = A c (t)[k + A(t)} + eA 2 c (t), 



t](t) 



I' 



A c (r)dT. 



(20) 



(21) 



(22) 



Eq. ([T5]) expresses the conservation law of energy for the 
electron tunnel ionized at the time t' . Eq. ([T9]) guar- 
antees that the electron returns to its parent ion as well 
as constraining the intermediate momentum of the elec- 
tron. Finally, (Eq. [20]) gives the energy conservation of 
the electron at the time t, when upon its return recom- 
bines with the core and releases a high frequency photon 

n. 

The terms A(t'), rj(t) and X(t) in the Eqs. (fTS]) . ([T9]) . 
and ([20]) , respectively, emerge from the nonhomogeneous 
character of the laser field and they vanish for the ho- 
mogeneous case, i.e when e = 0. For the homogeneous 
case, the solutions of the saddle equations are generally 
complex since Eq. ([15]) admits no real solutions, unless 
I p —> 0. This is a consequence of the fact that tunnel- 
ing has no classical counterpart. For the inhomogeneous 
case, however, it is not very upfront to constrain the 
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FIG. 1: Schematic representation of the laser electric 
field E(t) and the corresponding potential field A(£) for a 
monochromatic field denned by E(t) = Eq sin(cjt)e x , with 4 
cycles. The arrow indicates the classical times around which 
the electrons leave to the continuum and return at the field 
crossing, approximately. The pairs of the orbits are indicated 
by the labels Pair n , where n range from 1 to 8. Panels, 
a,b,c,d,e and f represent the first, second, third, fourth, fifth 
and sixth shortest pairs, respectively. The fields are normal- 
ized to A(t)/Ao and ~E(t)/Eo , where Ao and Eq are the field 
amplitudes. 
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III. RESULTS 
A. Quantum orbits 

In this section, we investigate the role of individual tra- 
jectories to the HHG cutoff for the nonhomogeneous case 
by performing a quantum-orbit analysis of the problem. 
The concept of the quantum-orbits is based on the fact 
that the solutions of the saddle-point equations can be 
related to the classical trajectories of the electron and, in 
addition, to obtain information on quantum aspects such 
as tunneling and interference. To get a better insight into 
the nonhomogeneous case, we employ a monochromatic 
field with E(t) = Eq s\n(ujt)e x , where e x is the polariza- 
tion vector along the x-axis. By using the relationship 
defined in Eq. (0 and applying some trigonometric iden- 
tities, the laser effective electric field ([T2]) and effective 
potential field (fT5]l along the electron trajectory read 

E tr (t) = £ sin(W)(l + 2esin(cjt)/cj 2 ), (24) 



A tr (t) = A cos(ut) + 2eA c (t), (25) 

respectively, where Aq — Eo /uj and 

A c (t) = Al s'm(ujt)/4uj - Alt/2. (26) 

In terms of the pondermotive energy of the homogeneous 
field, the drift momentum of Eq. ([23]) yields 

00 



k = -2y>%cos(ujt f ) + e( T f sin(ut') - 2U p t')(l =F Vs) 

(27) 
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FIG. 2: (Color online) Dependence of harmonic order on 
the release time t' and the recombination time t of the elec- 
tron for all given pairs in Fig. (pQ) and for the homogeneous 
field (e = 0). We consider an hydrogen atom, for which the 
ground-state energy is I p = 0.5 a.u., in a linearly polarized, 
monochromatic field of frequency uj — 0.057 a.u. and inten- 
sity / = 3 x 10 14 W/cm 2 . Panels (a) and (b) give the harmonic 
order as a function of the ionization and recombination times 
of the SFA model, respectively, while panels (c) and (d) depict 
the the harmonic order in terms of ionization and recombina- 
tion times of the classical calculations, respectively. The red 
dashed lines correspond to the harmonic cutoff. 



Based on the above equations, we solve the saddle 
point equations defined in Eqs. ([T8]) - ([20]) in terms of 
the ionization t' and recombination t times. For more 
close analysis, we restrict ourselves just to the solutions 
of the first 4 cycles of our defined monochromatic field, 
as shown in Fig. [TJ Classically, it is most probable that 
electron ionizes at the electric field maxima and returns 
to its parents ion at the electric field crossings. In FigHJ 
panels (a), (b), (c) and (d) depict the cases when the 
electron leaves at the field maxima and returns to the 
core at time about 7r, 27r, 37r, and 4tt later, respectively. 




0.75 0. 

Re[cot727r] 



FIG. 3: (Color online) SFA harmonic order as a function of 
the real part of the release time t' and the recombination time 
t of the electron for the same parameters as in Fig. [2] but for 
pairs l n and 2 m (where n — 1 — 6 and m = 1 — 5 ) of Fig. [T] 
Panels (a) (from right to left li — ► l^) and (b) (from left to 
right li — >• le) give the ionization and recombination times of 
pairs l n , respectively. Panels (c) (from right to left 2i — >> 2 5 ) 
and (d) (from left to right 2i — ► 2s) depict the ionization and 
recombination times of pairs 2 m , respectively. The dashed 
and solid lines correspond to the long and the short orbits, 
respectively. 

In Fig. [2j we plot the the harmonic order as function 
of the real parts of the ionization t' and recombination 
t times for the case with e = (panels a and b, respec- 
tively). In this figure, in comparison to the SFA model, 
we also present the classical solutions of t' and t (panels 
c and d, respectively). From Fig. [2j it is clear that the 
SFA resemble the classical calculations. Apart from that, 
both calculations show that the ionization and recombi- 
nation times corresponding to each cycle are identical. 

In Fig. [3j we isolate those solutions and plot the har- 
monic order in terms of the real parts of ionization and 
recombination times by considering the case when the 
electron leaves the atom at the times around tt/2 and 
returns to the core at the times nn (n = 2,3,4,5) (pan- 
els (a) and (b) for t' and t, respectively) and the case 
when the electron leaves the atom at the times around 
7r and returns to the core at nir (n = 3,4,5) (panels 
(c) and (d) for t' and t, respectively). For a given har- 
monic, there is always a shorter (solid line) and a longer 
(dashed line) travel time for the electron in the contin- 
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FIG. 4: (Color online) SFA harmonic order as a function of 
the imaginary part of the release time t' and the recombina- 
tion time t of the electron for the same parameters as in Fig. 
[2 but for pairs l n and 2 m (where n = 1 — 6 and m = 1 — 5 ) 
of Fig. [1] Panels (a) and (b) give the ionization and recombi- 
nation times of pairs l n , respectively, and panels (c) and (d) 
depict the ionization and recombination times of pairs Z m , re- 
spectively. The dashed and solid lines correspond to the long 
and the short orbits. 




"(b) 


SFA 










Re[cot'/27r] 



Re[©t/27c] 



FIG. 5: (Color online) Dependence of harmonic order on the 
release time t' and the recombination time t of the electron 
for all given pairs in Fig. [T] for a nonhomogeneous field 
with e = 0.003. We consider hydrogen atoms for which the 
ground-state energy is I p = 0.5 a.u. in a linearly polarized, 
monochromatic field of frequency uj — 0.057 a.u. and intensity 
I — 3x 10 14 W/cm 2 . Panels (a) and (b) give the ionization and 
recombination times of SFA model, respectively, while panels 
(c) and (d) depict the ionization and recombination times of 
the classical calculations, respectively. The red dashed lines 
correspond to the harmonic cutoff. 



uum, corresponding to the long and short trajectories of 
the pair. Such pairs of orbits coalesce at the maximally 
allowed harmonic energies, i.e at the cutoff. In here, the 
shortest orbit, i.e. the one that the electron leaves at the 
field maxima and returns at around time ir later, have 
the largest cutoff, at harmonic 47cj. It means that these 
pairs of orbits lead the cutoff of the HHG spectra, while 



the others pairs produce harmonics with lower energies. 

In Fig. HJ we present the imaginary parts of ion- 
ization and recombination times of the pairs shown in 
Fig. [3l Panels (a) and (b) depict the t f and t for nir 
(n = 2,3,4,5), respectively. Panels (c) and (d) demon- 
strate the t' and t for nir (n = 3,4,5), respectively. For 
the recombination times, Im[t] essentially vanishes be- 
tween the harmonic order for which the real parts Re[t] 
coalesce. Physically, this means that, in this region, the 
recombination is classically allowed. Beyond this region, 
Im[t] increases abruptly, which indicates that the clas- 
sically forbidden region has been reached. On the other 
hand, the imaginary part Im\t'\ of the start time of the 
electron is always non- vanishing. This is due to the fact 
that the tunneling has no classical counterpart. These re- 
sults show that both the imaginary and real parts of the 
tunneling and recombination times for pairs l n ( where 
n = 1 — 6 ) and 2 m (m = 1 — 5 ) given in Fig. [I] are 
similar. 
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FIG. 6: (Color online) Dependence of harmonic order on the 
release time t' and the recombination time t of the electron 
for the same parameters as in Fig. but for pairs m (where 
n = 1 — 6) of Fig. [U Panels (a) and (c) (from left to right 
li — ► 6i) depict the real part of the release and recombination 
times, respectively, while panels (c) and (d) show the imagi- 
nary part of the ionization and recombination times of pairs 
m, respectively. The dashed and solid lines correspond to the 
long and the short orbits. 

Now we move to the nonhomogeneous case and con- 
sider e = 0.003. In Fig. [5j we demonstrate the harmonic 
order as function of the real parts of the ionization t' and 
recombination t times for such a case. Like above, the 
SFA (Panels (a) and (b) for t' and t, respectively) resem- 
bles the classical calculations (Panels (c) and (d) for t' 
and t, respectively). The general harmonic cutoff is ex- 
tended in comparison to the homogeneous case, but the 
trajectories do not follow the same symmetry as shown 
in Fig. El 

To closely examine the orbits, we plot the harmonic 
order of the shortest pairs, i.e pairs nir (n=l to 5), in 
terms of ionization and recombination times (Fig. [6]). 
Panels (a) and (c) represent the real and the imaginary 
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parts of t' : respectively and, panels (b) and (d) represent 
the real and the imaginary parts of t, respectively. Unlike 
the homogeneous case, these pairs do not lead to the 
same cutoff. For the pairs corresponding to the electron 
leaving at the field maxima (Fig. QJa)), the cutoff is at 
around harmonic 38cj, while for the pairs corresponding 
to the electron leaving at the field minima the cutoff is 
at around harmonic 60co>. These results come from the 
fact that, for a given harmonic, the electron may tunnel 
with two possible momenta given by Eq. (|27|) . It appears 
that the electron has larger momenta if it tunnels from 
minima of the field and smaller momenta if it tunnels 
from maxima of the field. 




FIG. 7: (Color online) SFA harmonic order as a function of 
the real part of the release time t' and the recombination time 
t of the electron for the same parameters as in Fig. [5j but for 
pairs l n and 2 m (where n — 1 — 6 and m = 1 — 5 ) of Fig. 
[T] Panels (a) and (b) (from left to right li — > le) give the 
ionization and recombination times of pairs in 5 respectively, 
and panels (c) and (d) (from left to right 2i — »• 2s) depict the 
ionization and recombination times of pairs 2 m , respectively. 
The dashed and solid lines correspond to the long and the 
short orbits, while the pairs with dot dashed lines do not 
have the well-known shorts and long pairs. 

Figures and [8] demonstrate the dependence of har- 
monic order on the real and imaginary parts of ioniza- 
tion and recombination times, respectively, by consider- 
ing the case when the electron leaves the atom at the 
times around tt/2 and returns to the core at the times 
nix (n = 2, 3, 4, 5) (panels (a) and (b) for t' and t, respec- 
tively) and the case when the electron leaves the atom 
at the times around tt and returns to the core at nix 
(n = 3,4, 5) (panels (c)and (d) for t' and t, respectively). 

For the dominant pairs the cutoff become larger as we 
move from shorter pairs to the longer pairs. For instance, 
pair li, which is associated to the electron which leave 
at the times around 7r/2 and recombined at the times 
27r, has cutoff at harmonic 38a;, while the cutoff of pair 
16, which is associated to the electron which leave at the 
times around 7r/2 and recombined at the times 57r, is at 
harmonic 80co>. In contrast to the homogeneous case, in 
here, the pairs which associated with the electron which 



recombined after few cycle later from the time of its ion- 
ization lead to the larger cutoff. Furthermore, for some 
of the pairs (shown with dot dashed lines) like pair 1q 
and pair 23 the conventional concept of the short orbit, 
in which the electron leaves a bit later and returns a bit 
earlier, and the long orbit, in which the electron leaves 
a bit earlier and returns a bit later, does not have any 
meaning. In fact, for these pairs, if the electron leaves a 
bit early then return a bit early and if it leaves a bit later 
then returns a bit later. 
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FIG. 8: (Color online) SFA harmonic order as a function of 
the imaginary part of the release time t' and the recombi- 
nation time t of the electron for the same parameters as in 
Fig. but for pairs l n and 2 m (where n — 1,2,4,6 and 
m = 1, 3, 5) of Fig. [T] Panels (a) and (b) give the ionization 
and recombination times of pairs in 5 respectively, and pan- 
els (c) and (d) depict the ionization and recombination times 
of pairs 2 m , respectively. The dashed and solid lines corre- 
spond to the long and the short orbits, while the pairs with 
dot dashed lines do not have the well-known shorts and long 
pairs. 

In addition, the ionization time of the electron has even 
more interesting behavior as shown in Fig. [9] For pairs 
l n and pairs 2 rn (where n — 1 — 6 and m = 1 — 5 ), the 
ionization time t' moves toward tt as n and m become 
larger, until both solutions collapsed on each other as it 
demonstrated in Fig. [U In this figure, panel (a) depict 
the SFA calculation while panel (b) represent the classical 
calculations. For the SFA model, the ionization time 
as function of n and m moves more slowly toward tt in 
comparison to the classical calculations. In the classical 
calculations, the collapse at t' = tt is associated to the 
pairs I7 and 2q. In the SFA model, however, this collapse 
it is not exactly centered at around t' = 7r, instead it 
manifests itself by collapsing to its previous pair. On 
the other hand, both SFA and classical models give the 
same cutoff. It means the SFA will give a reliable HHG 
spectra while its yields will be affected quantitatively, 
which would not be a problem since SFA has the same 
limitation even for the homogeneous fields. 

In Fig. [lOl we demonstrate the harmonic order as func- 
tion of the real parts of the ionization t' and recombina- 
tion t times for e = 0.005. Panels (a) and (c) depict the 
t' and panels (b) and (d) show t, for SFA and classical 
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models, respectively. The SFA calculations are in good 
agreement with the classical model. In here, the general 
cutoff extended to the larger harmonic. In fact for the 
shortest pairs, the cutoff is at harmonic 70uj and for the 
longest allowed pairs the cutoff extends to harmonic 92cj. 
In comparison to the case with e = 0.003, for the shortest 
pairs (l n with n = 6), the cutoff from maxima and the 
minima of the field shift towards lower (harmonic 32co>) 
and higher harmonic (harmonic 70cj) as demonstrated in 
Fig. HJ 



100 




Re[cotV27t] 

FIG. 9: (Color online) Harmonic order as a function of the 
release time t' of the electron for the same parameters given 
in Fig. [5] but for pairs l n and 2 m (where n — 1,2,4,6 and 
m = 1,3,5 ) of Fig. [T] Panel (a) gives the real part of 
the release times based on SFA, while panel (b) depicts the 
classical ionization times of these pairs. The dashed and solid 
lines correspond to the long and the short orbits, while the 
pairs with dot dashed lines do not have the well-known short 
and long pairs. 

Furthermore, the collapse of the pairs from maximum 
and minimum of a given cycle happens even faster as 
shown in Fig. [121 As we discussed above for e = 0.003, 
the collapse occurs at Pairs I7 and 26, while in here it is 
at Pairs I4 and 23. It means that, as field becomes more 
inhomogeneous, i.e larger values of e, the electron does 
not return to the core if it follows the longer trajectories. 

We verified, without showing it here, that for larger 
values of e (like e = 0.02 ) even the longer trajectory of 
the shortest pairs does not lead to return of the electron 
to the core. Despite the fact that our model does not 
accommodate such large values of e, it suggests that for 
strong inhomogeneous field just the shortest trajectories 
will lead to recombination process. In fact, our TDSE 
calculations [22|, which based on the actual nonhomo- 



geneous field generated in the confined region of bow-tie 
nanostructures, show that the shortest trajectories rather 
than longest are contributing to the HHG spectra. 




1 2 3 1 2 3 4 



Re[cot727i] Re[cot/27i] 

FIG. 10: (Color online) Dependence of harmonic order on 
the release time t' and the recombination time t of the elec- 
tron for all given pairs in Fig. [TJ for a nonhomogeneous field 
with e = 0.005. We consider hydrogen atoms for which the 
ground-state energy is I p = 0.5 a.u. in a linearly polarized, 
monochromatic field of frequency uj — 0.057 a.u. and intensity 
/ = 3x 10 14 W/cm 2 . Panels (a) and (b) give the ionization and 
recombination times of SFA model, respectively, while panels 
(c) and (d) depict the ionization and recombination times of 
the classical calculations, respectively. The red dashed lines 
corresponds to the harmonic cutoff. 



B. Spectra 

In this section, we compute HHG spectra with Eq. (pp) 
and using the saddle point method developed in Sec- 
tion II. B. Fig. [13] presents HHG spectra for the case 
with e = 0.003. The spectra with yellow and blue colors 
represent the contributions from the shortest pairs with 
smaller cutoff, i.e. pairs 11,13 and I5, and the shortest 
pairs with largest cutoff, i.e. pairs I2, 14 and i6, while the 
red color demonstrates the contributions from all other 
orbits given in Fig. [TJ For all three cases, the HHG 
cutoffs are in good agreement with the trajectories rep- 
resented the previous section. The latter case, which 
corresponds to the longest pairs, leads to the largest cut- 
off. This is expected since for these pairs the electron 
has more time to accelerate in the field and return to the 
core with higher energy. On the other hand, for these 
pairs, the wave packet spreads too much in the contin- 
uum. Therefore, they lead to harmonic with lower yields. 
For the former case, which correspond to the shortest 
pairs, the harmonic yields will be large for both the cases, 
i.e. the one has largest and smallest cutoffs. The spectra 
with black color shows the total contributions, i.e from 
all the pairs given in Fig. [TJ 

We now zoom into the spectra of Fig. [13] to examine 
more closely the harmonic generated from these sets of 
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FIG. 11: (Color online) Dependence of harmonic order on the 
release time t' and the recombination time t of the electron 
for the same parameters as in Fig. [10] but for pairs l n (where 
n = 1 — 6) of Fig. [1] Panels (a) and (b) (from left to right 
li — >■ le) depict the real part of the release and recombination 
times, respectively, while panels (c) and (d) show the imagi- 
nary part of the ionization and recombination times of pairs 
l n , respectively. The dashed and solid lines correspond to the 
long and the short orbits, respectively. 
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FIG. 12: (Color online) Harmonic order as a function of the 
release time t' and the recombination time t of the electron for 
the same parameters as in Fig. [5j but for pairs li, I2, 14, 2i 
and 23 of Fig. [T] Panels (a) and (b) give the real part of the 
release and recombination times of pairs, respectively, while 
panel (c) show their imaginary parts. Panel (d) depict the 
classical ionization times of these pairs. The dashed and solid 
lines correspond to the long and the short orbits, while the 
pairs with dot dashed lines do not have the well-known short 
and long pairs. 



pairs (Fig. [14]). Up to harmonic 40cj, the shortest pairs 
with lower cutoff (li,l3 and I5) give the shape of the 
total spectra and from harmonic 40cj to 64cj the short- 
est pairs with higher cutoff (I2, 14 and le) dominate the 
shape of the total spectra. For higher harmonic the rest 
of the pairs lead the shape of the total spectra. 

Furthermore, it seems that both odd and even harmon- 



FIG. 13: (Color online) High-order harmonic spectra for 
hydrogen atoms (I p = 0.5 a.u.) and interacting with a 
monochromatic field of frequency uj = 0.057 a.u. and in- 
tensity / = 3 x 10 14 W/cm 2 for the case with e = 0.003. Yel- 
low and blue colored spectra show the contributions from the 
shortest pairs with smaller cutoff, i.e. pairs li, I3 and I5, and 
the shortest pairs with largest cutoff, i.e. pairs I2, 14 and 1q, 
while the red color demonstrates the contributions from all 
other orbits. Black colored spectra shows the total contribu- 
tions from all pairs given in Fig. [1] For clarity, all the HHG 
spectra are scaled. 
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FIG. 14: (Color online) Zooming into the HHG spectra given 
in Fig. M 



ics are present in the total spectra. At the lower regime 
of the spectra, one sets of harmonics are a bit more dom- 
inant than the others while at higher regime both odd 
and even harmonic have the same weight. 

In Fig. [15] we compute the HHG spectra for homo- 
geneous field (red colored spectra) and nonhomogeneous 
fields with e = 0.003 (blue colored spectra) and e = 0.005 
(black colored spectra), respectively. For all three cases, 
the HHG cutoff is in good agreement with the trajecto- 
ries analysis represented in the previous section and with 
the full 3D numerical calculations of Ref [21]. For the 
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case with e = 0, the cutoff is at around harmonic 45cj 
and for e = 0.003 the cutoff is at around harmonic 73uj. 
The field with e = 0.005 leads to the largest cutoff at 
around harmonic 93co>. For the latter case, there is less 
interference in the region corresponding the harmonics 
90u to 100cj. This behavior is expected since there are 
just two trajectories which contribute to such cutoff as 
shown in Fig. [TUl 




10 20 



40 60 80 100 
Harmonic order 



FIG. 15: (Color online) High-order harmonic spectra for 
hydrogen atoms (I p = 0.5 a.u.) and interacting with a 
monochromatic field of frequency uj = 0.057 a.u. and inten- 
sity / = 3 x 10 14 W/cm 2 . Red, blue and black colors depict 
the cases with e = 0, e = 0.003 and e = 0.005, respectively. 
All the HHG spectra are scaled for clarity. 



to a higher cutoff and the other to a lower one. Fur- 
thermore, we demonstrate that for an electron tunneling 
at the field maxima the tunneling time becomes larger 
when it recombines at a later time, while for the electron 
tunneling at the field minima the tunnel time become 
smaller as it recombines at the later time. In fact, as we 
go from the shorter pairs to the longer ones, the trajecto- 
ries from the minimum and maximum of the cycle moves 
towards each other until they collapse on each other at 
the field crossing. The nonhomogeneity character of the 
field determines this collapse. For larger nonhomogene- 
ity strength the electron only returns to the core if it 
follows the shortest pairs of trajectories. In addition, for 
some of the trajectories we can not define the conven- 
tional pair of long and short orbits. The former orbit 
corresponds to the case when electron leaves a bit ear- 
lier and returns a bit later and the latter case gives the 
trajectories in which an electron leaves a bit later and 
returns a bit earlier to the ionic core. Indeed, in here, 
if the electron leaves earlier it will then returns earlier 
and if it leaves later then it will return later. We also 
demonstrate that in the case of linear nonhomogeneous 
fields, both odd and even harmonics are present in the 
HHG spectra. At the lower harmonics one is a bit more 
dominant than the others while at higher harmonic order 
they both have equal weight. Within our model, we show 
that the HHG cutoff extends to the larger harmonics as 
a function of the inhomogeneity strength. 
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